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ABSTRACT 

This paper addresses the nonlinear interactions of a triad of initially 
linear stability waves. The triad consisted of a single two-dimensional mode 
at a given frequency and two oblique modes with equal and opposite spanwise 
wave numbers. The oblique waves were at half the frequency and streamwise 
wave number of the two-dimensional mode. Attention was focused on the 
boundary- layer transition at low frequencies and high Reynolds numbers. A 
five-zoned structure and low-frequency scaling were used to derive the 
nonlinear-interaction equations. In this study, we analyzed the initial 
nonlinear development of the waves; the results indicated that the two- 
dimensional wave behaves according to linear theory. Nonlinear interactions 
caused exponential-of-an-exponential growth of the oblique modes. This 
resonant amplification of the subharmonic depended on the initial amplitude of 
the two-dimensional wave and on the initial phase angle between the two- 
dimensional wave and the oblique waves. The resonant growth of the oblique 
modes was more pronounced at lower frequencies than at higher frequencies. 

The results presented herein are in good agreement with experimental results 
and offer new explanations of the observed process. 

1. INTRODUCTION 

Linear stability theory indicates that, for unstable boundary layers, a 
Tollmein-Schlichting wave at sufficiently small amplitude would grow and, 
subsequently decay along the boundary layer. This, by itself, cannot explain 
the transition to the turbulent state and the three-dimensionality of the flow. 



In the experiment of Klebanoff, Tidstrora, and Sargent (1962), a predominantly 
two-dimensional disturbance of prescribed frequency was introduced into the 
boundary layer by way of a vibrating ribbon. At relatively large Tollmein- 
Schlichting amplitudes, the experiment showed that a spanwise periodic three- 
dimensional structure evolved. This structure is characterized by spanwise 
alternating peaks and valleys, that is, regions of enhanced and reduced 
amplitude. The structure grows in intensity with downstream distance. The 
growth of the wave at peak positions is much larger than the original Tollmein- 
Schlichting growth and is in contrast to the slow Tollmein-Schlichting linear 
growth rate on the viscous scale. 

Linear stability theory cannot account for such three-dimensional waves 
growing more rapidly than the two-dimensional component of the initial 
disturbance. To explain this discrepancy, Craik (1971), among others, proposed 
a resonant-triad interaction. The triad comprises a two-dimensional wave and 
two oblique waves propagating at equal and opposite angles to the flow 
direction, such that the three waves have the same phase velocity in the 
downstream direction. Craik's (1971) temporal stability analysis of the triad 
showed that the nonlinear interactions can cause a large energy transfer from 
the mean shear flow to the oblique waves. This energy transfer would cause the 
oblique waves to grow rapidly and may be responsible for the rapid development 
of three-dimensionality in unstable boundary layers. Craik (1971, 1978, 1985) 
and Usher and Craik (1975) supported this analysis by approximating the 
interaction coefficients for the kind of interplay that might take place. 

Craik accounted for the resonance in one of the experimental conditions of 
Klebanoff et al. (1962), but not for that observed at other conditions. Since 
peak-valley splitting is not associated with the subharmonic waves, however, 
his results were criticized. 
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The presence of the subharmonic of the fundamental wave at the laminar- 
turbulent transition in a boundary layer is now generally accepted, having 
been well established by such experimental observations as those of Kachanov, 
Koslov, and Levchenko (1977), Thomas and Saric (1981), Saric, Carter, and 
Reynolds (1981), Kachanov and Levchenko (1984), and Corke and Mangano (1989). 
The appearance of the subharmonic causes the flow to become three-dimensional. 
This manifests itself in the formation of X-shaped structures arranged in 
staggered, alternating rows. Thus, aligned X-structure corresponds to the 
peak-valley splitting, whereas the alternating one corresponds to subharmonic 
resonance. In the transition described by Kachanov and Levchenko (1984), the 
parametric resonant excitation of a pair of oblique subharmonic waves by the 
plane fundamental wave is the main cause of the appearance of the low- 
frequency fluctuation. The oscillations with the subharmonic frequency reach 
considerable amplitudes only in the region of the upper branch of the neutral 
stability curve and behind it. Kachanov and Levchenko's (1984) observations 
indicated that three-dimensionality can come into play initially as a resonant- 
triad phenomenon; this revived interest in Craik's model. 

Attributing the occurrence of three-dimensionality in Klebanoff et al.'s 
(1962) experiment to spanwise differential amplification of Tollmein- 
Schlichting waves is supported by the theoretical work of Landahl (1972). 
However, the peak-valley splitting and the three-dimensionality of the 
disturbances may also be explained in the light of Craik's triad interactions 
by considering the two oblique waves to be at the excitation frequency but 
interacting with a two-dimensional wave at twice the excitation frequency. The 
peaks are thus given by the amplitudes of four waves: three waves at the same 

excitation frequency (one two-dimensional and two obliques) and a weaker two- 
dimensional wave at twice the frequency. The valleys are then given by the 
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amplitudes of two plane waves, one at the excitation frequency and one at twice 
the frequency. This explanation is supported by Klebanoff et al.'s (1962) and 
Cornelius' (1985) observation that the frequency tends to double at the 
valleys. With proper definition of fundamental and subharmonic frequencies, 
Craik's (1971) triad can apply to the transitions described by both Klebanoff 
et al. and Kachanov and Levenchko (1984). 

Of interest here are nonlinear interactions that arise from the continued 
growth of a resonant triad of initially linear instability waves. The triad 
consists of a single two-dimensional mode at a given frequency and two oblique 
modes with equal and opposite spanwise wave numbers. The oblique waves are at 
half the frequency and the streamwise wave number of the two-dimensional mode. 
Attention herein is focused on the technologically significant problem of 
boundary- layer transition at low frequency and high Reynolds numbers. Low- 
frequency scaling is used to derive the nonlinear equations. In the present 
analysis, nonlinearity first comes into effect in the common critical layer of 
these modes. This is consistent with the data of Klebanoff et al . (1962), with 
the numerical simulations of Croswell (1985) and Singer et al. (1987), and with 
Craik's (1971) finding that the nonlinear energy transfer takes place mainly in 
the vicinity of the critical layer. The amplitude equations are obtained by 
equating the velocity jump calculated from outside the critical layer to the 
velocity jump calculated from the solution within the critical layer. 

Two regimes of fundamental-subharmonic interactions can be identified in 
Kachanov and Levchenko's (1984) experiment. In the initial development of the 
three-dimensionality, the subharmonic is smaller than the fundamental, and the 
latter behaves as if the former is not present. Thus, although the subharmonic 
resonates, the fundamental is still given by its linear growth rate. In a 
second stage, where the subharmonic exceeds the fundamental, a back-reaction 
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process is observed in which the two waves behave nonlinearly. The initially 
negligible back-reaction on the two-dimensional wave allows the oblique modes 
to exhibit faster-than-exponential growth, while the two-dimensional mode 
continues to grow at its initial linear growth rate. This initial nonlinear 
interaction of the waves is the subject of the present study. 

The overall plan of this paper is as follows: the five-zoned structure 

and the scaling are discussed in section 2; the linear solution is obtained in 
section 3 (the nonlinear critical-layer interaction gradually evolves from a 
resonant triad of strictly linear, small growth rate and low-frequency 
solutions); in section 4, the critical-layer flow is considered; then, the 
amplitude equations are obtained in section 5, where the relevant expansions 
are worked out, and the solution that matches the solution outside the critical 
layer is then found by taking the Fourier transforms with respect to the cross- 
stream coordinate (the resulting solution is restricted to the case where the 
order of the oblique waves is not much larger than that of the two-dimensional 
wave); the solution of the amplitude equations in the nonlinear region is 
presented in section 6; the composite solution that matches the upstream linear 
solution is formed in section 7 (this leads to a different type of amplitude 
equation that involves upstream history effects and is different from that of 
Craik (1971)); results and discussions are presented in section 8; and a 
detailed comparison to available data is given in section 9. The experimental 
comparisons seem very encouraging, yielding both qualitative and quantitative 
agreement. The results presented herein explain the observed subharmonic 
resonance and the three-dimensionality of the flow. New features of the 
subharmonic resonance are discussed. 
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2 . FIVE-ZONED STRUCTURE AND SCALING 

The upstream flow starts as a resonant triad (fig. 1) of linear, spatially 

growing instability waves: a single two-dimensional mode of normalized 

frequency u and wave number a; and two subharmonic oblique modes of 

frequency u/2 , streamwise wave number nearly equal to a/2, and spanwise wave 

number ±|3. All velocities are normalized by the upstream velocity U^, 

2 

lengths by S* , time by SVU^, and pressure by pU^, where p is the fluid 
density and 6* is the boundary- layer thickness defined as 6* = y'yx/U^ . 

The flow under consideration here is that of a laminar boundary layer of 
Blasius profile (fig. 2). The development of the waves outside the nonlinear 
region follows the weakly nonparallel linear theory. In the nonlinear region, 
the unsteady flow is not affected by the boundary layer growth over the region 
in which the nonlinear interaction takes place. The experimental observations 
of Kachanov and Levchenko (1984), Saric et al . (1984), and Corke and Mangano 
(1987, 1989) have indicated that subharmonic resonance takes place in the 
vicinity and beyond the upper branch of the neutral stability curve. Goldstein 
and Durbin (1986) showed that upper-branch scaling applies over almost all of 
the unstable region, only breaking upstream in the close proximity of the lower 
branch. In this range of Reynolds number, the disturbance solution is a five- 
zoned structure (see Drazin and Reid 1981, Graebel 1966, Eagles 1969, Fraenkel 
1969, DeVillers 1975), as shown in figure 3, Zone I is the closest to the 
wall, where the unsteady viscous sublayer reduces to zero the effect of the 
streamline displacement at the wall. Zone II is an inviscid rotational zone of 
adjustment within which the critical layer III is induced. Zone IV is above 
zone II; it is an inviscid rotational region comprising most of the boundary 
layer, and it is the place where the disturbance nearly provokes a small 
quasi-steady displacement of the streamlines of the basic boundary- layer 
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flow. Above zone IV, the displacement also provokes the quasi-steady zone V, 
lying outside the boundary layer, in which the flow properties are of an 
inviscid irrotational type because the basic flow is almost the free stream 
U = 1 , V = 0 , where U is the velocity component in the flow direction and V 
is the perpendicular component in the y-direction. 

The normalized complex wave number a is small, and its imaginary part 
is smaller than its real part (Reid, 1965; and Goldstein, Durbin, and Leib, 
1987). Then each of the three modes has a critical layer at nearly the same 
transverse position y c , where the real part of their nearly common phase 
velocity c is equal to the streamwise velocity U. 

The mean boundary- layer velocity is given by the Blasius velocity Ug, 

where 

\ 2 

U B = Xy - 2 T 47 y + • • • , as y -» 0 . ( 2 . 1 ) 

The wall is located at y = 0 , and the constant X denotes the scaled Blasius 
skin friction (X = 0.332 in the present coordinates). 

We introduce the asymptotically small parameter o, where 

R = a 10 R . ( 2 . 2 a) 

This small parameter o can be related to the normalized frequency F*, as in 
Goldstein and Durbin (1986), by 

F* = a 12 ( 2 . 2 b) 

IT 

00 

Here, R is the Reynolds number based on the local momentum thickness, R is an 
order-one scaled Reynolds number and is a real quantity for time-periodic, 
spatially growing disturbances, and 

c » -§— (2.3) 

a 
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where c and a are order-one real constants. Goldstein et al. (1987) showed 
that the most rapidly growing modes correspond to the scaling 


Re <x = oa + O(a^) , 


(2.4) 


c o - cC 0 + °( a ) > 


(2.5) 


3 = a|3 , 


( 2 . 6 ) 


and 


y c = cY c + °(° ) - 


(2.7) 


where |3 is an order-one real spanwise wave number, Y is an order-one real 

c 

constant, and c Q is the phase velocity of the two-dimensional wave. The 
scaling of a and R given by equations (2.4) and (2.2a) is consistent with 
Reid's (1965) equation (3.128). 

3. LINEAR SOLUTION 

Outside the critical layer the unsteady flow is governed by the linear 
dynamics, as first pointed out by Haberman (1972), and the velocity field is 
given by 

u = U B (y) + Re |cA q (x 1 ) (y,x lt a)e lX + 8A(x ][ ) [u + (y .Xj ,a)e lZ 


U_(y,x 1 ,a)e lZ Je lX/2 j , (3.1) 


_ a_ . f . . iX . iZ -iZ, A iX/2l 
v = Ke l eaA o <t> o e + 8yA(e + e )<l>e I , 


and 


w = 6 Re A^W + e* Z + We ^ Z ^e^ 


iX/2 


where 


_ 4 

Xj = a x , 


(3.2) 


(3.3) 


(3.4) 
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(3.5) 


X E aa(x - act ) , 


Z E aftz , 

and 


(3.6) 


Y 



(3.7) 


The amplitudes A 0 and A are for the two-dimensional and the oblique waves, 
respectively. The symbols e and 5 are the measure of the amplitudes of 
the two-dimensional and the oblique instabi 1 i ty waves , respectively. Since 
the initial development of the instabi li ty waves is linear, we can take Afx^) 
to initially be a real quantity, but we allow the two-dimensional amplitude 
Ao( x i) to b e complex. Later on, because of the nonlinear effects, both A(x^) 
and A^xj) become complex. In equations (3.1) and (3.2), U ± and ff £ are 
related to $ through the relations 


and 



£-W = M> , 

Y ± 


(3.8) 



(3.9) 


D = i^ 


(3.10) 


By using the Squire (1933) transformation, the linearized equations for 


4> o and $ are 


and 




- i«R[(U - c q )(D 2 - <x 2 )* o 



(3.11) 


9 



(3.12) 


(D 2 - y 2 ) ^ = i yR[(U - c)(D 2 - Y 2 )^ - U”*] , 

where the primes denote differentiation with respect to the relevant arguments. 
The solutions of equations (3.11) and (3.12) are subject to the wall boundary 


conditions 


and 


$ = <f> = 0 at y = 0 
o J 


(3.13) 


$'=$^=■0 at y = 0 


(3.14) 


The complex wave numbers a and y and the phase velocities c and c are 

o 

given by 


a = aa + a 


. A' 

4 o 


iA ’ 
o 


(3.15) 


- a 4 a A' 

Y " °Y + 2l 


ac 


3 A [ ’ 


ia 


o 


and 


c = 


a c 


3 a > 

1 + 2 2 - L. 


la 


Equations (3.8) and (3.9) can be manipulated to produce 

3 


U cos e 4 W sin 0 = D4> - ^ 


U' . A' . 2 a 
+ -rj sin 0 , 


U - c iA 


where 


(3.16) 

(3.17) 


(3.18) 


(3.19) 


9 - sin 


in-1 §L 


(3.20) 


is the propagation angle of the oblique waves. 
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3.1. Solution in Viscous Sublayer (Zone I) 

For the viscous effects to come into play in this region, the scaling for 
y should be 

y - a 4 y , (3.1.1) 

where y is order one. In the classical theory, the wall layer thickness is 
order (aRc) -1/3 which is order 0(a 4 ), consistent with equation (3.1.1). 
Substituting (2.1) into (3.11) and (3.12), we obtain, up to the required order 
of approximation, 

[B 4 - iaR(Xa 3 y - c)D 2 ]$ = 0 , (3.1.2) 

^ ^ ~ - 
where represents or 4»; a represents a or y; c represents c q or 

c; and D is the derivative with respect to y. We expand $ as 

* = g Q + <7g x + O g 2 + O g 3 + . . • , (3.1.3) 

and with a and c expanded according to (3.15) to (3.18), we obtain 

( 8 4 + iRa C 0 D 2 )g 0 - 0 . (3.1.4) 

The solution of (3.1.4) not only must decay as y - ® but it also is subject 
to the conditions 

g o (0) - g o (0) = 0 . (3.1.5) 

The first-order solution is given by 

g o " m 2 " 1 + ^ ) ’ (3.1.6) 

m 2 - -iaRc Q , (3.1.7) 

where b Q is an arbitrary constant to be determined by matching with the 
solution at zone II. A similar solution is obtained for the oblique waves. 
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The corresponding velocity components are given by 
u - Re [cA o b o (l - e my )e‘ X + SAb cos 0 2 cos Z^l - e - ** y ^e i 


X/2 


+ . 




(3.1.8) 


Re 


-ia5 a if f my + e my - l)eA e lX 


- ia 5 y ^ ^My + e My - 1^6A 2 cos Z e iX/2 + 


(3.1.9) 


and 


w = Re j^2ib sin Z sin 9^1 - Me My )sAe iX/2 + . . .j , (3.1.10) 


where 


M = -iyRc . (3.1.11) 

3.2. The Tollmein Region (Zone II) 

As pointed out by Goldstein, Durbin, and Leib (1987), the solution for 
this region is obtained by introducing the scaled transverse coordinate 


(3.2.1) 


directly into equations (3.11) and (3.12) and by using equations (2.1) and 
(3.15) to (3.18). The solution is obtained in the form 


<t> 


= a [( c 1 + ac 2 + ° 2ci i) + ( c 2 + oC 4 + ct2c 6 ) Y 1 + ct4f • (3.2.2) 


5 >] 


As will be shown, the matching requirements give c^ = C 2 - C3 =* 0; therefore, 
the solution can be written as 


$ = o(X + oa)Y + <7 F(Y,<|>) 


(3.2.3) 
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where a is an order-one constant that depends on ct, and F satisfies the 
relation 


where 


0 - * rrV;) - 1 x2y(y c * Y > ■ 


2 2 
X 

C 


(3.2.4) 


(3.2.5) 


and c = XY 


Since equation (3.2.4) is singular at Y - Y c , F can be discontinuous 
across Y c , and F* denotes the solution above or under this point. 

Integrating equation (3.2.4) gives 

(■*<¥,♦) - F(0) * M c || V 2 t Y C [(Y - Y c )(ln|Y - Y c | + i**) 

:(>" Y c * *♦■)]} - 4TJT y3 ( Y c + \ ¥ ) ' < 3 - 2 ' 6 > 


+ Y 


where the normal velocity is continuous across the critical layer. The 
constants of integration and <j>* are, in general, a complex function of 

. The value of F as Y - 0 is F(0), which is determined via matching 
with the solution in zone I. The corresponding velocity components are given 
by 


u 


= Re eA Q [(X + oa) + a 3 F(Y,<j>)je* X + 26A cos Z^(X + aa)cos 9 + [j““ 


x sin 0 tan 0 o(X + aa)Y + crcos 9 F'(Y,<j>) 

4 U’ . ~ . - viv . J iX/2 

-i- a u - sin 9 tan 0 F(Y,<|>)Je , 


(3.2 .7) 


Re | -eA Q iaa 2 ^(X + oa Q )Y + o 3 F(Y,<j>)J e*" 


- SAiya 2 [(X + ca)Y + o 3 F(Y,(|))]2 cos Z e lX/2 j , (3.2.8) 
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and 


w = jzi SA sin Z sin 9^(X + aa) - g ^ c a(X + aa)Y 

+ * 3 F(Y,<|>) - a 4 F(Y,$)]e iX/2 j. (3.2.9) 

3.3. Solution in Zones IV and V 

Miles' (1962) solution, which is uniformly valid for y = 0(1) and y » 1 
in the limit as y,a - 0, is given by 


where 


D$ U' 


+ 0(a 5 ) 


$ u - c (u - era 


1 ~2 

-^r~2 + + afl 1 + a 0 2 + . . . , 


a(l - c) 


(3.3.1) 


(3.3.2) 


0 = - 
0 


(i - c y 


A O A n 

HL-- c) 2 (i - c) 2 


A r\ A n 

L(i - c ) 2 (U - c ) 2 J 


dy , 


(3.3.3) 


a 


1 (i - c ) 2 J 


(U - c) Q Q dy , 


(3.3.4) 


and 




(u - cr 


20 


(1 - c) 


1 A 2 

+ Q 

2 o 


dy . 


(3.3.5) 


By substituting equations (3.3.1) to (3.3.4) into the classical "inviscid 
function" (Lin, 1955, p. 37), which is defined as 


I = 


cD$ 

U’$ - (U - c)D$ 


(3.3.6) 


and by inserting U = Ug and equations (3.15) to (3.18) into the result, 
expanding for small o, and then using equation (2.1), we obtain 
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where 

„,+ t - ~ \ U'c , „ -\-2 ocX Tt „ ~r 2 - 2 r 12 3^1 

W (y;a,c,a) - -x— (1 - ac) - _ 4 1 J t + 2 ocJ 2 + a c J 3 + ^ y + 0 (y )J 

a (1 — ac) 

+ o 2 ac\ J 4 + acJ 5 - ^3 y + 0(y 2 ) + a 3 c -§ In y - |g ^ + °( y ) • 

4X . .X X 

(3.3.9) 

Coefficients to J g are defined in the appendix. 

3.4. Matching Solutions 

Matching the solution in zone I and the solution in zone II is easily done 
by matching the velocities. The inner solution for the two-dimensional wave 
(eqs. ( 3 . 1 . 8 ) to ( 3 . 1 . 10 )) as y * ®, reduces to 

u = $0 b Q e lX , (3.4.1) 

v = Re [-iao 5 b o (y + . . .]e lX . (3.4.2) 


For the solution in zone II (eqs. (3.2.7) to (3.2.9)), as Y - 0 , we have 


u ” ke £x + oa + ct 2 F'(o) + . 

. .].« . 

(3.4.3) 

v - && ^-ia 5 a Q F'( 0 ) +■ . . 

F . 

(3.4.4) 

Matching the u-velocities gives 



b o " X . 


(3.4.5) 


whereas matching the v-velocities gives the imaginary part of F( 0 ) as 



(3.4.6) 
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Matching the velocities of the obliques waves produces similar expressions, 
with c replacing c 0 and y replacing a. 

Matching the solution in zone II and the outer solution in zones IV and V 
is most easily achieved by using the inviscid function. The inviscid function 
for the solution in zone II is obtained by substituting equation (3.2.3) into 
equation (3.3.6) and re-expanding: 


2 U 3 
" - X * ° 




Y - Y 


In 


— " - ** - W 1 ~ l XY c Y < 2Y c + Y > 


+ 0(cr 4 ) (3.4.7) 


for Y > Y c , where A$ denotes either - 4> + ) or ($ Q - <j>g^. Matching 
equations (3.3.7) and (3.3.8) shows that 


A' 


a -f- cf 2 2 9m-p- + 0(o 4 ) , 


(3.4.8) 


X §-( 1 - ac)“ 2 - 1 - 
a 


+ a 


acX 


(1 - ac) 4 


Ol + 2a8 h + ° z ~ z 1 3 ) 


5 (r >" Y c 

7 5 2 
48 .2 


aj 

F_(0) XY 2 
- XY + 4X a 

\ 

X 

a 

0 ) 

c 


+ a 2 acX^J^ + ocjg^ 


= 0 , 


(cos 0 + 


A' 


- 2 - 
Y V. 


cos 0/ A 


-2 

c y 

Y^(A<D) + 


-.2 c 
cX 


^2cRy 


and 


A o 


- 2 - 
a u 


T Y_(A4> ) + 

9 C rj 

2c \ 

o 


-2 

a 


2 c 0 Y c X y2c o Ra 


cXY c Y 2cR Y 


t \ 1 

*1 Uf 

A 


O A 


(3.4.9) 


(3.4.10) 


(3.4.11) 


'initial 

Equations (3.4.8) and (3.4.9) are dispersion relations which determine a 
and y in terms of c (or in terms of the scaled Strouhal number ac). Since 
their coefficients are all real, they are consistent with our original 
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assertion that a and c are real quantities. In fact, it follows from 
equations (3.7) and (3.4.8) that 

3 - ^ a + — <?9m J 9 + 0(o 4 ) . (3.4.12) 

1 y/3 o 

This shows that 3 and a satisfy the usual long-wave length, smal 1 -growth- 
rate resonance condition to within the order of the detuning. To a first 
approximation, equation (3.4.9) shows that c 0 and a satisfy the usual 
long-wave length, small-growth-rate dispersion relation 

5 „-J- < 3 - 4 - 13 > 

The imaginary part of the dispersion relations (i.e., the matching between 
zones II and III) produces the real parts of equations (3.4.10) and (3.4.11). 
For the initial linear growth, A<j> is real. The amplitude A is real, and A q 
is initially complex, but its imaginary part remains equal to its initial 
value. Nonlinearity causes both amplitudes to become complex. Since the 
amplitudes appear as order a in the dispersion relations, the imaginary 
parts of A' /A or A^/A q are balanced by the imaginary part of the 
corresponding A<{>. Therefore, equations (3.4.10) and (3.4.11) in their complex 
form are valid for both the linear and nonlinear growth regimes; they relate 
the (slow) growth rates of the instability waves A' /A and A^/A o *° 
phase jumps A<)> and across the critical layer. To determine these 

latter quantities, we must consider the flow in the critical layer. 

Equations (3.2.3), (3.2.6) to (3.2.9), (3.4.5), and (3.4.6) show that the 
flow field in zone II can be written as 
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— — <J> sin 2 0 

XY - c 


(cos Z)6Ae 


v = -o 2 XY Re [iaeA Q e iX + 2iy(cos Z)6Ae lX/2 ] + . . . , (3.4.15) 

w = -2S(sin 0 sin Z) Re — — z iAe* X ^ 2 + . . . , (3.4.16) 

XY - c 

and 

p = acX Re ^eA Q e iX + 2 (cos 0 cos Z)SAe iX/2 ] + ... , (3.4.17) 

where 

f . Mc (Y ♦ Y c In | Y - Y c l >Y c )~f Y 2 ^ * ^ . (3.4.18) 

This solution clearly becomes singular at the critical level Y = Y c and has 
to be rescaled in this region. 

4. CRITICAL LEVEL 

As indicated in the preceding section, the solution obtained for zone II 
(eqs. (3.4.14) to (3.4.17)) becomes singular at the critical layer. The 
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appropriate transverse coordinate in this region, which matches the linear 
growth rate, is 


Y - Y_ 


(4.1) 


Reid (1965) predicted the critical layer thickness of the order (aR) -1 ^, which 
is O(ct^), in agreement with the scaling in equation (4.1). Introducing 
equation (4.1) into equations (3.4.14) to (3.4.17) and expanding the result for 
small n shows that the critical-layer solution should expand as follows: 


3, - 6 + v - 8 1 f - \2 . -2 f -1 c c r 2 

u-c = aXr) + a p„Y^n + o p c 2 (n) + 5a u 0 + 6c u_, + 5u n + 6au 1 + 6a u 


c c 


-1 


, 3 (o) (o) 2 (o) 3 (o) 

+ 5a u^ + eUp + ecu^ ' + Co u| } + ecr u^ ’ + 

v = -ea -2 aXY Re iA e*^ - Sa -2 YXY 2 cos Z Re iAe*^ 2 + 

p n 1 r* 


- 2-1 23 

w = 5a w_ 2 + Sa w_ 1 + Sw Q + SoWj +■ 5a vi^ + 5a w g + 


and 


p - eacX Re A e lX + 5acX cos 0 2 cos Z Re Ae 1 * + 
o 


where we have set 


and 


4- 

v = a v 


2 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 


3 2 


(4.7) 


The linear analysis indicates that Uq, u!j\ and u° are independent of n and 
that 

U 2 " U 2 = W 0 = W 2 " 0 * (4>8) 
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The linear solution ((3.4.14) to (3.4.18)) shows that the functions u° and 
U 3 are discontinuous and singular across the critical layer. The other 
functions u_ 2 , u_i, u^, w_ 2 , w_ lf and wj are singular at n = 0 , and 
therefore, the linear solution is invalid. The full nonlinear momentum 
equations are thus needed to obtain a valid solution at the critical level. 

The nonlinear terms in the critical-layer solution will balance the 
velocity jump at the same order as in the linear solution when 


c 3 < 7^ 8 . (4.9) 

By substituting the expansions (4.2) to (4.5) into the full momentum 
equations, the latter can be expressed in terms of the scaled variables 
xj, X, Z, and n as 


Du = - <xP Y + a 3 P , a 8 P- , 3P y 

~ \ a Xj n Lj 

and the continuity equation can be written as 


(4.10) 


au v + v- + w„ + a u, = 0 , 
A r| Zj 


(4.11) 


where we have set 


3 . - 8 _ . _3„ _3_ . n... 3_ _3 U £ 

'l w “ 9n‘ 


D = S < U - c > 3X + 5 ^ * 0 u 3J7 * B * 3Z - 0 H .-2 


(4.12) 


and 


u = (u , v , w ) 


(4.13) 


The Haberman parameter 



l 


R 


(4.14) 


is order one. 
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It is also convenient to work with the equation for Z-component vorticity 
o), which can be written as 

I _2 g g_ \ 

aDu - pl a u^w- + atow^ - a av^w^ - a v^w x J . (4.15) 

The first stage of the development of the waves, which occurs when the 
amplitude of the oblique waves does not exceed the order of that of the two- 
dimensional wave, is considered here. In this case all modal amplitudes are 
of the same order so that 

6 = e (4.16) 

Substituting (4.16) into (4.2) to (4.5), combining terms of equal power, and 
substituting the resulting expansion into (4.10) and (4.13), gives 

<g o u_ 2 = (r 5 ) p 2 2 cos Z Re iAe lX/2 , (4.17) 


U 1 ’ -™0x, - V 1 X - “ P 3X + V 


^cV U -2X + ^ U -2x, 


\ V 2 + 
XY u 

C 


E , 


&0 u 3 - - 5u 2Xl - V 3 x - pSX + P 2 x J - lx + i^J u < 


(4.18) 




-2 


2 + u -2 J ax XY c E a - + w -2^ az 


, -2 x 

\ n ) 

2 + U -2 ’ 


(4.19) 


J?_ w = ^-|-P 2 sin Z Re Ae lX/2 

0 -2 y l 


W 1 - -“Vc^ W -2X - X ^-2 Xl - ^3Z - CW 0x 1 ’ 


(4.20) 

(4.21) 
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o -2n 


<£ u 1 - = -cu n - + X(Sw.™ 
=^o In OnXj \ 1Z 


&o u 3n ■ -tfcfaoiX - S 'oz) - X 


in 


Onx, 


; * xy c*-2R E 


w -2 u -2x) * a,1 c 2 "-2X ' 

(4.22) 

W -2Z ’ 

(4.23) 

" V^X * f ' u -2Xn) ’ 

(4.24) 

5, 3Z - V c( u -2nn 1 0 & E 


fz2 ^ 



and 


aU -2X + ** W -2Z " 0 ’ 

“ U 1X + % + * W 1Z + u -2 Xl " 0 ’ 

« U 3X + % + ^ W 3Z + u o Xl " 0 ’ 


(4.26) 

(4.27) 


(4.28) 


where 

~ aX ^ lx “ H 7^2 

on 

and 


(4.29) 


E - i(A Q ae lX + 2 Ay cos Z e lX/2 ) . (4.30) 

These equations ultimately determine the amplitude equations. Equations for 
velocity components at other o orders are not displayed here, since they do 
not influence the amplitude equations. The y-momentum equations indicate 
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that the pressure derivatives with respect to n are zero up to the 
order . 

5. AMPLITUDE EQUATIONS 

The set of equations (4.17) to (4.28) must be solved subject to the 
transverse boundary conditions so that they match the outer solution, equations 
(4.2) to (4.5). Equations (4.17) and (4.20) indicate that the lowest order 
solutions are given by 

u_ 2 = 2 tan 0 cos Z fie £iQ(n .x^e*^ 2 ] (5.1) 

and 

w_ 2 = 2 sin Z fie Qffj.XjJe*^ 2 . (5.2) 

where Q satisfies 


= yXc(sin 0 cos 0)A , 


subject to the transfer boundary condition 


We have set 


Q - 0 as n - ±«> . 


E n = Xf > - H bl for ” " *< 2 ■ 


(5.3) 


(5.4) 


(5.5) 


3n 

Equations (4.22), (4.23), (4.25), and (4.26) are manipulated to obtain 

"in ■ ■' u oSx 1 - at ‘c Y c(r-2) iiX + - Xfiu -2x in - X, -2x 1 sin 6 

and 

Sto U 3S ■ “^(f" ^ u inX - D \) * x Vc cos 9 E - - Dt, o) 


(5.6) 


- X (' iu SnX 1 - D+W 3 - ¥ c fie Eu -2Ss) 

- (*-2 D+U -2 - u -2 D+W -l)- - ^ U 2nx. ' 15 - 7) 

n i 
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where 


u* - u n cos 0 + w r sin 0 for n - 1,2,3 


and 


(5.8) 


D + = (cos 0)0 - (sin 0)a fy . 


8X 


(5.9) 


The functions u*~ and u*- must be expressible in the form 


_3 

> 

m,n=0 


S C5,x 1 )e I[(n/2 J x+(m - 1)zl (5.10) 


and 


u 3n ‘ & L Cx (n- 1 )e ,[(n/2,Xt( ”- 1,Z1 • 
m,n=0 


.CD 


where Q* ' satisfies the relation 

1 9 1 


(5.11) 



( 1 ) 

1,1 



sin 0 , 


Q 


(3) 

1,1 


satisfies the relation 


(5.12) 


E X!1 - « r n 2 J ♦ w c [« sin 6 *0 4 * V A cos 9 ] ' < 513 > 

13 1 

and Q satisfies the relation 

E 2 Q 2^0 ' cos e [ iS tan20 ^ " XY c*( A tan 6 ' 2iv *c A o cos 6 )] ■ (5,14) 

Here, the asterisks denote the complex conjugate. 

Examination of equations (3.4.14), (3.4.17), and (4.2) to (4.5) indicates 
that only ut- is discontinuous across the critical layer. Therefore, it 

dp 

follows from equations (3.4.10), (3.4.11), and (5.11) that these solutions 
must satisfy 
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cX 2 i 

-2 

Y 


X 2 A 


cos 6 + —^- 5 ) A' - 
, cos QJ ^ Rj 


(5.15) 


and 


n (3 ) .- 2cX Z i 

Q 2 , 0 dn • ~-2 


a 


A' - - 
0 2V2 R 


X 2 A „ /A'Y 

*- - “> 1 


cos 9 


(5.16) 


in order to match the discontinuous 0(a^^) terms in equation (3.4.14). 
Equations (5.15) and (5.16) arise from the requirement that the changes in the 
fundamental and subharmonic components of the velocity across the critical 
layer as calculated from the external solution be the same as when they are 
calculated from the internal solution. This latter solution is obtained by 
solving equations (5.3) and (5.12) to (5.14) and by substituting the result in 
equations (5.15) and (5.16). 

To obtain the inner solution, we introduce the following normalized 
variables: 

x = x 1 - x Q , (5.17) 


n 

c 


2 H 


- -3 ’ 
aXc J 


(5.18) 

(5.19) 


where x c is the origin of the nonlinear region. Introducing these into 
equations (5.3) and (5.12) to (5.14) and taking the Fourier transforms with 
respect to n gives 


Q - 2 ir(sin 0)A(x)H(-k)e 


hk 3 /3 


(5.20) 


jsn( 1 ) 
9 Q 1,1 
3x 


ir ein 2 fi dA 

— sin 9 JT 

ac ax i 


hk 3 /3 


f H(k)dk , 
-k 


(5.21) 
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q( 3 ) 

^1,1 


4irY c hk 3 /3 
— — e 


ok 


§Aj ^-) A (x)A*(x) 
c / ° 


-2hk 3 /3 „ 

, A. 


k^HfkjJdkj + i y A(x)H(-k) 


(5.22) 


5(3) 

'2 ,0 


hk 3 /6 
= e 

, _ \3 

±r&_ [_§_] a 2 


p 00 p 


Xac 3 v Y / 




V 

00 ^ 

-00 J. 


( h /3)(k^ k |)_[( hk 2) /e ] 


2ir(3 2 Y 


C .2 


x 6(k x + k 2 - k 3 )dk x dk 2 dk 3 - — 3-e A 


ac 3 y 


W) 


/6 


k i ^ -k l^ e dk 


V 

+ — liru H(-k) 
yc C 


(5.23) 


where the notation 


oo 

F(k) = J e _ikn F( n )d n 


(5.24) 


is used to denote the Fourier transform of any function F(n) so that its 
inverse is given by 


00 

F (n) - ^ J e ikn F(k)dk 


(5.25) 


Since the integral over n is just the zero wave number component of the 
Fourier transform, then equation (5.15) gives 

o / \ 2 


c 


cos 9 * ssh) 1 A * V + irA * 4 o ' t 5 - 26 ) 


where 


*Y_ , - N 2 


\ ^ ( t -) 


(5.27) 


and 
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r - » ^ ^Lr sin 2 e . 
a X 


(5.28) 


Solving for Q 2 ^ shows that the nonlinear interaction term makes no 
contribution to the velocity jump integral of equation (5.16), so it follows 
that the fundamental amplitude A q is still determined by the linear growth 
relation 



which can be integrated to give 


ff_R? 

8 X 3 + 


2y/2 R 


where 



v 


(5.29) 


(5.30) 


it R 2 


8 ,3 


X" 2\/2 R 


(5.31) 


This shows that the fundamental is determined by the linear growth relation 


(k +ik. )x 
o 1 o i' 


a o e 


(5.32) 


where a° is the complex amplitude of the fundamental instability wave at 
x = 0, and kj is the initial wave number detuning factor, which can be chosen 
as an initial condition. This shows that there is no back-reaction of the 
oblique mode on the two-dimensional mode, and that, consequently, the two- 
dimensional wave follows the linear theory. 

With the fundamental wave given by equation (5.32), the oblique waves' 
amplitude equation (5.26) is written as 

i - 1 k o A * Di v* ■ (5 - 33) 
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where 


D - 


3_ 

10 



(5.34) 


and the following relations are used to the order of approximation of the 
analysis : 


c - f - XY . 4 

X c a 


a - y. cos y - C = XY„ > -r . (5.35) 

Equation (5.33) can be split into magnitude and phase equations, respectively: 


dlAJ 

dx 


I k | A| 

5 o' 1 


^-= D 
dx 


D 


k x 
o 


e sin(t oi + k.x - 2^)|A| 


k x 


cos(v|/ q + k { x - 2q#) , 


(5.36) 

(5.37) 


where ij/ is the phase angle of A and is the phase angle of a°. The 

solution of equations (5.36) and (5.37) is subject to the condition that it 
match with the upstream linear solution. 

6. DEVELOPMENT OF AMPLITUDES WITHIN THE NONLINEAR REGION 
In the nonlinear region, k Q , k. , and R can be taken as constants. The 
amplitude of the two-dimensional wave follows the simple linear growth relation 
of equation (5.32). The solution for the oblique waves' amplitude equations 
(5.36) and (5.37) depends on the initial detuning factor. In section (6.1), we 
present an analytical solution for the case of perfect tuning, that is, kj = 0. 
The effect of the detuning factor will be discussed in section (6.2.) 

6.1. Perfect Tuning 

For kj - 0, the system of equations (5.36) and (5.37) can be solved 
analytically to obtain 


28 


(6.1) 


|A| 


\ ) D 

O 

^ J K 

a o 


P + 2 ln 


1 + c 4 e 


- 4 (D/k o ) 


1 + c. 


where 


v 


P - e 


and the effective phase angle t defined as 

c 


(6.2) 


is given by 


where 


te , +o - 2* 


sin(i|/ ) = i 

' T e' r + 1 


r = c. e 
4 


- 4 (D/k o ) 


a° P 
o 


( 6 . 3 ) 


( 6 . 4 ) 


( 6 . 5 ) 


The constants a and c 4 are to be determined from matching with the outer 

linear solution x ■* -<*>, and c d can be related to the phase angle y at 

00 


x - 0 by 


1 + sin y 


eo 


'4 1 - sin 14/ 


( 6 . 6 ) 


eo 


The first exponential factor in equation (6.1) is the linear growth rate. 
The next exponential factor represents augmentation or suppression of the 
growth above the linear growth, depending on the initial phase angle teo 
(which determines C4). Note that C4 ranges from 0 to ®. Two values of teo 
are of particular interest. If f eo - ir/2 , 04=*® and equation (6.1) reduces 
to 
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(6.7) 


1*1 - a ex p(i k o ; ) e,<p (ir | a o| e ° ) ■ 
If i|r go - 3ir/2, c 4 - 0 and equation (6.1) reduces to 



For - 3ff/2, the second exponential factor in equation (6.8) represents an 
explosive, exponent ial-of-an-exponent ial growth. Thus the resonant-triad 
interaction is a very powerful mechanism for augmenting the growth rate of the 
oblique instability waves that propagate at an angle of 60° to the downstream 
direction. 

The amplitude of the oblique waves in the nonlinear region (eqs. (6.1) and 

(6.4)) at various initial phase angles q» eQ is shown in figure 4 for a° - 1 

and R = 0.8 R , where R is the upper-branch scaled Reynolds number, 
up up 

Figure 4(a) shows that the amplitude grows explosively for q< eo = 3ir/2, as 

given by equation (6.8), The trend is reversed when q> eQ =» ir/2 (eq. (6.7)). 

The phase of the amplitude is shown in figure 4(b). If t = ir/2 or 3ir/2, 

the effective phase angle q» e remains at its initial value. If q» eQ = 0 

or it, the effective phase angle develops rapidly to an asymptotic value of 

»|» = 3ir/2. 

T e 

6.2. Effect of Detuning Factor 

Since perfect tuning may not always be achieved, we discuss here the 
effect of the initial detuning factor lq. The magnitude of the plane wave's 
amplitude, equation (5.32), is independent of the detuning factor. The 
detuning factor causes the phase angle of the plane wave to vary linearly 
in x. The effect of the detuning factor on the oblique waves is obtained by 
solving equations (5.36) and (5.37) numerically for several values of kj; the 
results are shown in figure 5. Only large values of lq , compared to k 0 , can 
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influence the solution. Figure 5(a) shows that detuning only moves the 
location of the resonant amplification downstream; however, detuning causes the 
effective phase angle to vary in the upstream linear region, as figure 5(b) 
indicates. But the effective phase angle reaches its optimum value of 
sin i|» e = -1 some distance downstream, which allows for the explosive growth to 
occur. This can be understood by noting that the solution of equation (5.37) 
is independent of that of equation (5.36), and that the former can be written 
as 


where 


dx. 


dx 


k. 


2D 


k o* 

e cos x. 


(6.9) 


x i “ ^o + ^i x ” ^ ' (6.10) 

For large x values, equation (6.9) reduces to that of perfect tuning. Thus, 
for <|» eo = 270°, if kj = 0, sin = -1 is achieved immediately at x = 0. 

As kj increases, the optimum angle is achieved farther downstream. The 
nonlinear contribution to equation (5.36) is dependent on kj in sin(xj). 
Since sin(xj) reaches the optimum value at some downstream location that 
increases with kj , the only effect of kj on the amplitude is to move down 
the location of resonance amplification. Since for reasonable values of kj 
the effect of detuning is thus negligible, the solution at kj - 0 can be 
considered a good approximation. Therefore, in the following sections, where 
the composite solution and the results are presented, we assume k[ - 0 and 
use the analytical solution (eqs. (6.1) to (6.5)). 

7. COMPOSITE SOLUTION 

The overall instability-wave growth is strongly dependent on both the 
nonlinear effects discussed here and the weakly nonparallel effects. In order 
to incorporate both effects into a single formula, we use the composite 
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expansion as suggested by Van Dyke (1975). The composite expansion is formed 
from the present inner solution Aj n and the slowly varying outer solution 
A ou t- Here the composite solution for the case of perfect tuning, that is, 
kj =* 0, is presented. In section 7.1 the composite expansion in the unsealed 
coordinate x, which is the dimensional streamwise coordinate normalized by 
the local boundary- layer thickness, is obtained. In section 7.2 the integrals 
that depend on the boundary- layer variations are evaluated and the solution is 
presented in terms of the unsealed boundary- layer coordinates. 

7.1. Composite Expansion of the Two-Dimensional Wave and Oblique Waves 

Consider first the plane wave. The outer linear solution can be written 
as 


L k o< x 2> dx 


A . = A (x_)e 
o,out o v 2 ’ 


(7.1) 


where x 2 is the slowly varying boundary- layer variable. The origin of the 
disturbances is at x = 0 (see fig. 2). As x 2 -» x^, this outer solution can 


be approximated by 


o 


k x 

A . a° e 0 = A . . 

o,out o 0,1/0 


(7.2) 


where 


rX 

f k € 

Jn O 


a° = A°(o)e 0 ' (7.3) 

Consider the inner solution of the plane wave (eq. (5.32)) as x ♦ -®; the 
solution can be written as 

k x 

A. -+ a 0 e 0 = A . . , . (7.4) 

in o o,in/out 

The composite solution is given according to van Dyke (1975) as 
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(7.5) 


A A 

» o. in o,out 

A r ” A 

t>, in/out 

Thus, for the plane wave, the composite solution is given by 


ct 4 f k (x )dx 
Jn o v 2 


A = A (o)e J 0 
o,c o v ’ 


(7.6) 


and we have approximated A q (x 2 ) by 


A o (x 2 } = A o (x = 0) 


(7.7) 


The composite solution for the plane wave (eq. (7.6)) is thus the same as its 
outer solution. This is an expected result since the two-dimensional wave in 
the nonlinear region is still given by its linear solution. 

With respect to the oblique waves, the outer solution is given by 

X 

'o 


A . » A°(x„)e 
out v 2 


(4/5)<j 4 f k dx 
Jo O 


(7.8) 


As x_ -* x , this outer solution can be written as 

6 O 


A . + a e 
out 


(4/5 )k x 


in/out 


(7.9) 


where 


(4/5)a 4 f ° k dx 
Jo O 


a = A u (o)e u . (7.10) 

The inner solution for the oblique waves, equations (6.1) to (6.4), can be 
written as 


( 4/ 5 )k x _ 
A in - a e g(x) 


(7.11) 


where 
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0 


sin (v o - 2vj/ ) = 


r - 1 
r + 1 


and 


-4T 

r = c 4 e 


As x ■* -®, 

(4/5)k Q x 

A in * a e = A in/out * 


(7.14) 


(7.15) 


(7.16) 


The composite solution 

A. A . 

, in out 
A - T 

c ln/out 

is then given by 

(4/5 )<7 4 J k Q dx 

A = A°(o)e 0 g(x)e 1 '*' 

c 

where A°(x 2 ) has been approximated by 

A°(x 2 ) = A°(x - 0) . 


(7.17) 


(7.18) 


(7.19) 


We allow the coordinate system origin shift x q to occur naturally by 
rewriting equation (7.13) as 



pX 


T = 

D 

•S 

V 



d _o> 


V 


ax . 


(7.20) 


Substituting for 



from equation (7.3), we obtain 
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A"(x = 0) 


a 4 f k dx 

J 0 0 

D e U dx 


(7.21) 


The constant c 4 , equation (7.15), is given by 


1 + sin . 

T oi 


1 - sin if 


01 


(7.22) 


where i|» Q . is the initial phase angle of the fundamental. 

7.2. Composite Solution in Unsealed Boundary-Layer Coordinates 
Here the solution of the amplitude functions in the unsealed boundary- 
layer coordinates is presented. The composite solutions (eqs. (7.16) and 
(7.18)) can be integrated once the boundary- layer variation with x is given. 
Let 


I 



(7.23) 


where k (x_) is given by equation (5.31) and x is normalized by the local 

0 M 

boundary- layer thickness $*. If 6* is the boundary thickness at x = 0, we 
can write 


* 



(7.24) 


where R. is the scaled Reynolds number at x = 0. If x is the distance 
1 " 

between the leading edge of the flat plate and the origin of the disturbances 
x = 0, then according to Schlichting (1955), we can write 



(7.25) 
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where x^ is the dimensional streamwise coordinate. Equation (7.23) can be 
integrated to obtain 


I = ct 4 R. 


24 x 3 




x^ + 1 (x^ + 4) 


- 4 


X 2 1 


4 -s/2 R. 


ln(x T + 1) + 


x f + 1. 


(7.26) 


where 


As in Goldstein and Durbin (1986), 
by the dimensionless frequency F* , 

F* = 


the scaling parameter 
defined as 
wv 12 



(7.27) 

a can be replaced 

(7.28) 


and equation (7.26) can be written as 


I = - 


R?F* 2 
ir i 

24 x 3 


(x* + 4) tfj + 1 - 4 + — — ln[x^ + 1) + , x ^ 

J 4n/2 v/F 7 L x T + 1. 


The composite solutions are now written as 


A 


A 0 

o| 


0 


I 

e , 


| A I - A 0 e( 4/5 > T ; 


g(T) 


(7.29) 

(7.30) 

(7.31) 


and 


Sin i); 


-4T , 

c 4 e - 1 

e -4T , 

c 4 e +1 


(7.32) 


where g(T) is defined in equation (7.12). From equation (7.21), T can now be 
written in terms of the boundary- layer coordinates as 
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T = (F*) 


0.333 


> 

o 

R. 

1 0| 

l 


V 


cX' 


1 ± DeW . 


n fi3/2 
(1 + x‘ ) 


(7.33) 


The coordinate can be replaced by the local Reynolds number since 

\2 


xt *fe - 1 


(7.34) 


The scaled initial amplitudes A° and A° can be related to the unsealed 

amplitudes A° and A° „ by 
r n o,n 

A° = F* 10/12 A° (7.35) 


and 


A o _ f *10/12 a o 
o,n o’ 

8. RESULTS AND DISCUSSION 


(7.36) 


Results are presented here for the development of the fundamental and the 
subharmonic waves as given by equations (7.30) to (7.32). The solution of 
these equations is dependent on the initial levels of the fundamental and 
subharmonics and on the initial phase-difference angle between the two- 
dimensional mode and the oblique modes. These factors, among others, will be 
discussed here. 


8.1. Effect of Initial Level of the Fundamental 
Figure 6 shows the amplitudes of the fundamental and the subharmonic waves 
normalized by the initial level of the fundamental. The normalized frequency 
F* = 10 -5 , and the initial phase angle is 3ir/2. The unsealed initial 

level of the subharmonic is fixed at 10 _ ®, whereas the unsealed initial level 
of the fundamental is varied between 10 - ^ and io - 2.5 < Since the fundamental 
behaves linearly, as indicated herein, the figure shows that the amplification 
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of the fundamental is independent of its initial level. The development of 
the fundamental is also independent of that of the subharmonic. Kachanov and 
Levchenko (1984) have indicated that the development of the fundamental wave 
does not depend on the subharmonic development in an initial region of the 
amplification of small priming subharmonic fluctuations. Their results 
indicate no back-reaction on the fundamental up to a fundamental's initial 
level of 0.218xl0 -2 , which is consistent with the present results. 

In this analysis, the subharmonic does not initially react back on the 
fundamental. This allows the fundamental to grow according to the linear 
theory. On the other hand, the subharmonic grows as the exponent ial-of-an- 
exponent ial , as figure 6 shows. This is consistent with the observations of 
Kachanov and Levchenko (1984), Saric, Kozlov, and Levchenko (1984), and Corke 
and Mangano (1989). This growth of the subharmonic is dependent on the growth 
of the fundamental, which acts as the catalyst for the resonant growth. The 
growth of the subharmonic waves increases with increasing fundamental initial 
amplitude, and as figure 6 indicates, the subharmonic saturates with the 
fundamental's decay. 

8.2. Effect of Initial Level of the Subharmonic 
An interesting feature of the present analysis is that the subharmonic 
amplitude is given by a linear function in A. The ratio A/A° and the 
amplification rates are thus independent of the initial level of the 
subharmonic. The role of the initial amplitude of the subharmonic is clearly 
shown in the experiments of Saric et al. (1984) and Kachanov and Levchenko 
(1984, fig. 14(c)) in the case of simultaneous ribbon excitation at Tollmein- 
Schlichting frequency and its subharmonic. In the experiment, the initial 
level of the subharmonic was varied sinusoidally. This resulted in an 
amplitude A of the amplified subharmonic, which varied sinusoidally as the 
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initial disturbance and was proportional to the initial amplitude of the 
phase-locked component of the ribbon excitation. This indicates that the 
growth of the subharmonic is linearly dependent on the initial level of the 
subharmonic, which is consistent with the present analysis. 

8.3. Effect of Initial Phase-Difference Angle 

The magnitude of the fundamental's amplitude, which is governed by a 
linear process, is not dependent on the phase angle, as equation (7.6) 
indicates. On the other hand, the amplitude of the subharmonic, 
equation (7.31), is dependent on the constant C 4 . This constant is 
determined via the initial phase angle, q^o i > by using equation (7.22). Thus, 
the nonlinear process is dependent on the initial phase angle y 0 i between 
the fundamental and the subharmonic. 

The effect of i|/ 0 [ on the growth of the subharmonic is shown in figure 7. 
The effective phase angle te controls the sign of the nonlinear terra in the 
subharmonic's amplitude equation (5.34). For y 0 i » ir/2 or 3ir/2 , 
equation (7.31) reduces to 

| A | - A 0 e< 4/5)l e ±T , (8.3.1) 

where the negative sign is for i)r . - ir/2 and the positive sign is for 

dr . = 3ir/2 . The nonlinear effects can thus amplify or reduce the growth rate, 
T oi 

depending on the initial phase angle. The amplitude of the subharmonic at 
R/R i - V2 is shown in figure 7(b) as a function of t oi • The figure shows 
that, except in the vicinity of t|/ ^ = ir/2, the amplitude is close to its 
maximum value, which occurs at = 3ir/2. 

If . 3ir/2 , equation (7.22) shows that c 4 = 0 and thus 

equation (7.33) gives sin * = -1 irrespective of x. If y Qi - ir/2, c 4 is 

infinity and is also a fixed value independent of x, equation (7.33) 

gives sin - 1. If is not close to ir/2, c 4 is finite and, as x 
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increases, >|» reaches an asymptotic value sin w = -1. Thus, except for 
6 6 

«l» . - ir/2, sin t tends to be -1 as x increases; these results are 
demonstrated in figure 8. The sign of the nonlinear term in the oblique waves' 
amplitude equation (5.36) is determined by sin t Q . . Thus, except for 

s ir/2, the amplitude of the subharmonic is close to its value of 

^oi = 3lT ^ 2 ’ as was ^own in figure 7(b). 

In free-shear flows, the dependence of the subharmonic amplification on 
the initial phase difference was observed and studied by several researchers. 
The analysis of Monkewitz (1988) and the numerical simulations of Patnaik, 
Sherman, and Corcos (1976) and Riley and Metcalfe (1980) for vortex pairing in 
two-dimensional mixing layers, showed that the pairing process depends on the 
phase difference between the fundamental and subharmonic instability waves. In 
the experiment of Zhang, Ho, and Monkewitz (1985) for a two-dimensional shear 
layer under bimodal excitation, significantly different merging patterns were 
observed as a result of changing the initial phase difference between the 
fundamental and the subharmonic. For a round jet under fundamental-subharmonic 
excitation, Mankbadi (1986) showed that the initial growth rate of the 
subharmonic is highest when the two waves are in phase. The effect of phase 
difference on the subharmonic's peak was found to depend on the Strouhal -number 
pair. For a round jet under two frequency excitations, Mankbadi, Raman, and 
Rice (1990) demonstrated that the interaction between the fundamental and 
subharmonic waves depends on an effective phase angle identical to the one 
defined here in equation (6.3). 

For boundary layers, Klebanoff and Tidstrom (1959) showed that the 
disturbance grew to a large amplitude, which appears not to be qualitatively 
reproducible. This can be explained by the present analysis, which indicates 
that these high peaks may or may not occur, depending on toi ■ Since no 
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control was exercised over t(/ 0 i in Klebanoff and Tidstrom's (1959) experiment, 
it could vary at random. The analysis herein showed that if Y 0 j was optimum, 
peaks were observed, and if y 0 i was close to ir/2, then peaks were not 
observed; thus, the variation of the initial phase-difference angle can explain 
why the observed peaks were not reproducible. 

For the case of resonant amplification of controlled subharmonic priming 
oscillations, Kachanov and Levchenko (1984, p. 231) observed that resonance 
occurs at a given phase shift Yr between the crests of the fundamental wave 
and the subharmonic. They also reported that Zelman and Maslennikova's (1982) 
numerical results showed the same feature. Zelman and Maslennikova (1982) 
calculated the growth of the waves at various initial phase shifts between the 
subharmonic and fundamental. The results of their calculations show that when 
an initial phase of the subharmonic minus the resonant phase, measured on the 
cycle of the subharmonic, is 90° (or 180° measured on the cycle of the 
fundamental), the subharmonic not only fails to grow, but it also damps with 
approximately the same rate as the amplification at y r . Other things being 
equal, the damped subharmonic lags behind the subharmonic having y 0 i = by 
an order of magnitude or more in amplitude. These numerical results are 
consistent with those obtained herein (fig. 7(a)) for the cases of Y 0 i - ir/2 
and 3ir/2 , where the difference between the two values is ir on the 
fundamental cycle or tr/2 on the subharmonic cycle, as in the numerical 
simulations. 

8.4. Nonlinear Function 

The nonlinear function g in equation (7.32) is the ratio of the 
nonlinear amplitude to the linear amplitude. This function, defined in 
equation (7.12), can be written for tj/ 0 j - 3ir/2 as 
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where is the initial Reynolds number. Thus, 2. < 0 corresponds to the 

location of the origin of the disturbances (e.g., vibrating wire) ahead of the 

upper branch of the linear neutral curve, and 2 > 1 corresponds to a location 

of the origin of disturbances beyond the upper branch of the neutral curve 

(i.e., in the linear stable region). Other parameters being the same, 

equation (8.1) indicates the nonlinear interaction increases with the initial 

Reynolds number, R. . Also, as shown in figure 6, the nonlinear interaction 

increases with A 0 , which reaches its maximum at the upper branch. These two 
o 

factors explain why the observed subharmonic resonance occurs at or beyond the 
upper branch. 

The effect of frequency on the nonlinear function is shown in figure 9, 
which indicates that the nonlinear effects increase considerably with 
decreasing the frequency. Note that the argument of the exponent I in 
equation (8.1) also increases with decreasing frequency, as given by 
equation (7.29), Thus, the nonlinear interactions are quite sensitive to the 
normalized frequency F* . Since low values of F*, correspond to higher 
transitional Reynolds numbers, the present nonlinear phenomenon is particularly 
significant at high-Reynolds-number flows. 
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8.5. Nonlinear Upper Branch of Neutral Curve 
Because of the nonlinear effects, the subharmonic can continue to grow 
beyond the upper branch of the neutral curve. The amplification factor for the 
oblique waves, based on the linear theory, is 0.8k o , whereas 0.81 is that 
which accounts for the variation of the boundary- layer thickness. However, the 
total amplification rate of the oblique modes is given by the sum of the linear 
and nonlinear amplifications. For q/ 0 | =. 3ir/2, the nonlinear upper branch of 
the neutral curve for the oblique modes can be obtained from equation (7.31) by 
setting 

0.81 + T - 0 , (8.3) 

which can be written as 


A 0 R 4 


0| Up ^ 


1 + 0 . 5x' 

(1 ♦ x+) 3/2 


The solution of this equation determines the location of the nonlinear upper 
branch. This equation indicates that the nonlinear neutral curve is dependent 
on the initial Reynolds number as well as on the initial level of the 


fundamental. The first term in equation (8.4) is the linear growth rate of the 
subharmonic, which includes the small nonparallel flow effects. The second 


term is the nonlinear contribution to the neutral curve. As A° -> 0. the 

o 

nonlinear contribution vanishes, and the nonlinear curve coincides with the 


linear one. 


Since the nonlinear contribution in the upper-branch equation (8.4) is 
always positive for nonzero values of A°, nonlinear effects extend the 
critical Reynolds number. Thus, depending on A° and R q , the upper branch of 
the neutral curve can be eliminated. Goldstein and Durbin (1986) have pointed 
out that nonlinear critical layers eliminate the upper branch of spatially 


43 



growing Tol Imein-Schl ichting waves. Figure 10 shows both the linear and the 
nonlinear upper branches of the neutral curves for the subharmonic. The 
boundary- layer growth is accounted for in both curves (A° = 0.001 and 

R./R^p = 0.8^. The figure shows that at high frequencies the two curves 
approach each other. But at low frequencies the nonlinear upper branch's 
Reynolds number is much larger than that of the linear one. The upper-branch 
Reynolds number can be increased further by increasing the initial level of 
the fundamental. Thus, at low frequencies, the upper branch is practically 
eliminated, so an initially linear stability wave can encounter nonlinear 
effects that cause a state of suspended, nondecaying waves, as in the fully 
turbulent case. 

9. DETAILED COMPARISON WITH EXPERIMENTS 


The results presented in section 8 are in qualitative agreement with the 
available data. Further detailed comparisons with the data of Kachanov and 
Levchenko (1984) and Corke and Mangano (1987) are presented in this section. 

In figures 11 to 13 the amplitudes of the two-dimensional and oblique waves 
are compared with the corresponding experimental data. In calculating these 
amplitudes, the initial amplitudes and Reynolds numbers were taken to be the 
same as in the experiment. 

Figure 11 shows the calculated amplitudes at F* = 137xlO - ® compared with 
Kachanov and Levchenko's (1984) experimental data. The figure shows excellent 
agreement in the initial region. Further downstream, for R > 600, the 
measured amplitude of the oblique waves is less than the calculated one. 
Figures 12 and 13 show the calculated amplitudes at other frequencies and 
initial conditions compared with the corresponding data of Corke and Mangano 
(1987). These figures also show an excellent agreement for the initial 
development of the waves. The plane wave grows almost at its linear growth 
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rate. But both theory and observations show considerable growth of the oblique 
waves beyond that of the linear theory. At points downstream, the measured 
amplitude of the oblique waves is less than the predicted one. Since the 
amplitude of the subharmonic exceeds that of the fundamental, the disturbed 
flow enters a higher amplitude regime, and rescaling of the amplitudes is 
required. The present theory thus successfully predicts the initial resonance 
of the oblique waves. When the amplitude of the oblique waves exceeds that of 
the two-dimensional wave, a second mechanism comes into effect that reduces the 
growth of the oblique waves. 

10. CONCLUDING REMARKS 

The nonlinear interactions between a two-dimensional fundamental stability 
wave and a pair of oblique subharmonic waves of equal and opposite propagation 
angles were studied by using low-frequency scaling. Attention was focused on 
the initial stage of the nonlinear development of the waves. The present 
analysis indicated that the two-dimensional fundamental mode growth was given 
by the linear theory, whereas the subharmonic oblique modes grew explosively as 
the exponent ial-of-an-exponential . The fundamental acted as a catalyst for the 
growth of the subharmonic; therefore, the nonlinear growth of the subharmonic 
increased with an increase in the initial level of the fundamental. 

The development of the subharmonic was found to be dependent on the 
initial phase-difference angle between the two-dimensional fundamental mode and 
the oblique subharmonic modes. For the case of perfect tuning, the optimum 
initial phase difference was 3ir/2, measured on the cycle of the fundamental. 

If the initial phase difference was changed it from its optimum value, the 
subharmonic not only failed to grow, but decayed at a rate almost equal to its 
amplification rate for the optimum initial phase-difference angle. 
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The nonlinear growth of the subharmonic increased as the nondimens ional 
frequency F* decreased. Thus, the present triad interactions are stronger at 
lower frequencies. Keeping other parameters the same, the nonlinear 
interaction increases with Reynold's number. Nonlinear effects tended to 
increase the upper-branch Reynolds number of the oblique modes. At low 
frequencies, the upper branch of the neutral stability curve for the oblique 
modes was practically eliminated. 

Detailed comparisons with available data indicated that the present 
analysis successfully predicts the initial resonance of the oblique waves. If 
the subharmonic was not much larger than that of the fundamental , both theory 
and observations indicated that the fundamental grew according to the linear 
theory while the subharmonic grew at rates much higher than that predicted by 
the linear theory. The strong agreement between theory and experiment here 
has established the origin of the subharmonic resonance and of the observed 
three-dimensionality of the flow at boundary- layer transition. 
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APPENDIX 


In this appendix we list the numerical coefficients that appear in 

(3.3.9). 
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Figure 1 . — Resonant triad. 
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Figure 6. — Amplification curves for two-dimenslonai funda- 
mental A 0 and subha/monic oblique waves A at different 
Initial amplitudes of fundamental wave |A £ n |. Dimension- 
less frequency F* - 1 0 location of origin of disturbances 
R i/Rup - 0.8; phase angle ¥<* - 3 tc/ 2; and subharmonic 
inital level, A ° * ^ 0 ~ 8 . 
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Figure Qr - Dependence of development of the phase 
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Figure 7— Effect of initial phase-difference angle V oi on 
amplification of subharmonic. Dimensionless frequency 
F* * 1 0 - ' 5 ; location of origin of disturbances R j /R up - 0 .8; 
Initial amplitude of fundamental | Ag.nl * 0.00 1 ; and Inital 
amplitude of subharmonic A£~ io~® 


Normalized 



Figure 9. — Nonlinear amplification of subharmonic g for 
various values of dimensionless frequency F\ Initial 
amplitude of fundamental |A§, n |* 0.001; initial amplitude 
of subharmonic Afj » 10** 8 ; phase anjgle Vqi * 3jc/ 2; and 
location of origin of disturbances R j/R up * 0.8. 
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Figure 10. — Upper branches of neutral stability curves for 
nonlinear and linear development of oblique waves. 

Initial amplitude of fundamental jA g in | * 0.001 ; phase 
angle * 3 tl/ 2; and location of origin of disturbances 
R j /R up ■ 0.8. 
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Flflure 12.— Comparison of calculated amplitudes of fun- 
damental and subharmonic streamwise velocities wth 
experimental data of Corke and Mangano (1987) at 
F*« 93x10“®. 
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Figure 1 1 .—Comparison of calculated amplitudes 
of fundamental and subharmonic streamwise 
velocities with experimental data of Kachanov 
and Levchenko (1984) at F 4 « 137x1 0 ~ 6 . 
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Figure 1 3.— Comparison of calculated amplitudes of fun- 
damental and subharmonic streamwise velocities with 
experimental data of Corke and Mangano (1 987) at 
F 4 - 82.7x10 “ e . 
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